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Abstract 

We consider the Nuclear Modification Ratios in heavy ion colhsions, Rcp and Kaa, in the re- 
gion of intermediate transverse momentum, and study the dependency upon the constituent quark 
composition of the observed hadron. Adopting a two component recombination/fragmentation model, 
validated by experimental information from STAR and PHENIX, we show that a clear distinction is 
predicted for the /o(980) between the assumptions of ss or diquark-antidiquark content. 
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1 Introduction 



Some of the most important results obtained so far at RHIC are related with the high-px spectrum in 
heavy ion collisions. A very strong suppression of high-p^ particles [T] is commonly interpreted in terms 
of energy loss of the energetic partons when traversing a very opaque medium (see e.g. [21 [3] for recent 
reviews); two particle correlations offer the first jet-like measurements in heavy ion collisions, with some 
surprising features which await explanation [U [5] . 

On the other hand, unexpectedly, the protons yield dominates over the pion yield at intermediate 
values of 1 GeV< p± ^ 6 GeV [6], in contrast with any other result in experiments involving hadrons. 
This increase in the relative yield of protons is in fact common to other baryons and the medium effects in 
this px-region have been found to be driven by the number of constituent quarks of the particle studied. 
In particular, the amount of suppression has been experimentally found to depend on quark counting 
rather than on particle mass. Similarly, the azimuthal momentum anisotropics, measured by the second 
order, V2, of a harmonic expansion in the azimuthal angle, were found to exhibit scaling with the particle 
quark number (for a recent review of the data see, e.g. [7]). An economical model to describe these 
features of intermediate-px particles is the recombination or coalescence model O |TOj . 

The overall picture is that one exploits the different behavior exhibited by parton fragmentation in the 
hard sector and by hadronization in the soft sector: the cost to produce a particle in the fragmentation of 
a high-p_L parton increases with the number of constituent quarks, while recombination gives essentially 
the same spectrum due to its additive nature in momentum. Due to the dominance of the recombination 
component, baryon relative yields are enhanced in central A + A collisions with respect to peripheral 
collisions or p + p collisions. Predictions for multicharmed hadron abundances in the coalescence regime 
have been recently presented in [TTj . 

Other models, based on a hydrodynamical description of a heavy ion collision describe reasonably the 
soft part of the spectrum (see e.g. [E]). Models based on baryon formation mediated by string-junction 
configurations can also account for the observed enhancement in the ratio baryon / meson [13] . These type 
of models have been applied before to SPS data to explain the strong enhancement of strange baryon 
production |141ll5j. Understanding the dynamical origin of these effects, which suggest a strong interplay 
between hadronization, thermalization and parton shower evolution, needs of additional observables. 

The validity of the coalescence model in a sizeable range of pj_ may serve a dual purpose. It confirms 
the formation of an unconfined state in central heavy ion collisions and it opens up the exciting possibility 
to use the unconfined state to probe the constituent quark composition of hadrons, thus discriminating 
between conflicting models. One proposal in this direction [16] was to measure the azimuthal anisotropy 
of the possible resonances which were interpreted as pentaquarks a couple of years ago. The difficulty is 
that, on one hand these resonances ~ if real resonances at all - are difficult to measure and, on the other 
hand, the azimuthal anisotropics is a very small, percent level, signal. 

In the present paper, we consider the nuclear modification ratios, Rcp and Raa, as experimental 
observables with large signals, and find that they are suitable to measure the valence quark content of 
the resonances and allow to clarify their nature. 

We take, in particular, the case of the /o(980), proposed to be a 4-quark meson in [T7] and already 
measured in heavy ion collisions [I8j. Being a well known resonance, abundantly produced and easy to 
identify, the /o(980) provides a unique opportunity for these types of studies. 

We consider the two alternatives; /o(980) = ss [19] or /o(980) = [iTj and show that they can 

be very well distinguished by data at RHIC. We consider the extension of our analysis at LHC and find 
rather similar results. Other interesting cases are the /o(1370) and /o(1710) for which interpretations as 
qq or tetraquark states have been advanced, in [20] and in |21j respectively. 
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According to the tetraquark interpretation, /o(1710) has the same composition as /o(980). The Rcp 
and Raa predictions for this state are in fact essentially the same as for /o(980). The /o(1370) is supposed 
to be a state [nn][nn] {n = u,d), a recurrence of the ct(500) meson in the multiquark scheme. A brief 
discussion of the expected signal for /o(1370) is given at the end of this paper. An open problem in this 
context is how a glueball hadron would scale, e.g the /o(1500) as proposed in [20l[21]- Also lacking are 
predictions for a molecular state, such as the /o(980) in [22] . 



2 The model 

The soft part of the spectrum is assumed to be formed by the decay of a deconfined state of constituent 
quarks which recombine to form the observed hadrons, following valence quark counting. The hard part 
is described by perturbative QCD with a simple implementation of the energy loss of the high energy 
partons traversing the medium. 

We summarize in this Section basic formulae and parameters used for the two components of the 
hadron production model, recombination and fragmentation, following Ref. [8]. Notice that gluon degrees 
of freedom are not included in the original model, they have been discussed in [9] but the results are 
basically unchanged and will not be considered here. 



Recombination. The basic ingredients are the constituent quark or antiquark distribution functions, 
which we write according to: 

=^^-la- Pa{Pa,T) • faipa, <^a) , (1) 

where ja is the fugacity of the constituent and fa is a function that takes into account the azimuthal 
dependence of the distribution for non-central collisions. For recombination from a non-relativistic fluid 
at temperature T, 50, would simply be EadV, dV being the volume element, and pa the Boltzmann 
distribution. For the case we consider, namely recombination from a relativistic, boost-invariant fluid 
[23| . 6^1 is proportional to the volume element of the 3-dimensional, boost-invariant hypersurface defined 
in 4-space by: 

= {t, X, y, z) 
t"^ = — = const.; t = rcoshr/, z = rsinhry 

= normal to hypersurface = (cosh?], 0_l, sinhr/) . (2) 

Accordingly: 

= Tdrj d^x±{u ■ pa) ■ (3) 

and x± = (x, y) is the 2-dimensional transverse space. 

The Boltzmann distribution is replaced by the Touschek |24j invariant phase-space distribution: 

Eg {vPa) . ^. 

e T ^ e T ^ (4) 
We denote by the local 4-velocity of the fluid and by v± the transverse expansion velocity: 

= (cosh r/ cosh ry_|_, sinhr/_|_ cos 0, sinh r/_|_ sin (/>, sinhrycoshry_|_) 
v± = tanh rj± . (5) 
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Denoting by Ua and 4>a the constituent rapidity and azimuth in the transverse plane and by nia 



+ the transverse mass, we have: 



{V -pa 



^a,± cosh rjj_ cosh(r/ — y) — p± sinh ri± cos( 



(6) 



In the coalescence component, hadrons with transverse momentum Pi are formed by recombination 
from the thermal spectrum of partons, which share the momentum of the hadron equally, in the most 
favorable situation. 

We specialize, for simplicity, to the case of a qq meson and use the notations of Ref. [8], hadron 
variables are denoted with capital letters. The meson momentum distribution at mid-rapidity is given 
by: 



Em 



dNM\ _ I "^VAf 



(- 



lalbC, 



M 



j rdrj r±dr±d(/) (u ■ P) 



(7) 



where Cm is a multiplicity factor and the hadron momentum. 

We first neglect the azimuthal dependence by setting fa = fb = ^- The integrals over rj and (p give 
rise to approximate exponentials and one finds, in conclusion [8j: 



dNM 



d'^P±dy 



tAj_ 

y=Q {2ttY 



P± sinh r]± 



k2{Pi 



(8) 



A± = tt/Oq is the transverse area of the partonic medium at freeze-out, r the hadronization time and M_l 
is the transverse mass of the hadron. We have introduced the Bessel functions: 



1 fn roc 

Iq{z) = - #6"^™'^'^; Ki{z)= (i7?coshr? 6"^^°'^'^'' 
71" Jo Jo 



(9) 



and the short-hand notation for a N-quark hadron: 



kN{Pi 



cosh r]± 



N 



a=l 



Pi 



For baryons and 4-quark states we find, similarly: 

dNB 



d^P±dy 



y=0 



CBM^^^2-^alblch 



■(27r)3 



4g 



d^P±dy 



y=0 



(2^ 



Pj_ sinh r]x_ 
T 

Pj_ sinh ?7_L 



HPl) 



(10) 



(11) 



(12) 



Following [8] the degeneracy factors are taken to account for the spin of the final hadron only and 
we do not distinguish between A and T,^, so that: Cm = 1> Cb = 2 for protons, Cb = 4 for A -|- S'^ and 



G 



4q 
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Until now we have considered central collisions. For peripheral collisions a further, common fugacity 
factor for all parton species 7periph) is included and the transverse size of the medium is scaled with the 
impact parameter, b, according to: 



Po Po ^2 (13) 



^A 

where Ra is the radius of the nucleus. 

We list in Table [1] the values of the parameters assumed for central and peripheral collisions at 
RHIC. We follow here mostly Ref. [8], where it is shown that these parameters reproduce reasonably the 
intermediate p± cross sections at RHIC. 

Table 1: Parameters of the hadron production model in Au+Au @ RHIC, Pb+Pb @LHC and p+p. Numbers in parentheses 
in the last column refer to peripheral collisions with & = 12 fm at RHIC and & = 14 fm at LHC. 





lu,d 




1s,s 


Tperiph 


tAt 
(fm^) 




eo 

(GeV-i/2) 


^coll 


Au+Au @RHIC 


1 


0.9 


0.8 


0.7 


1.27-10^ 


0.55 


0.82 


1146 (26) 


Pb+Pb @LHC 


1 


1 


1 


0.7 


11.5-10^ 


0.68 


2.5 


3643(82) 


p+p @RHIC(LHC) 


0.4 


0.4 


0.12 (0.4) 




13.2 (119) 


0.55 (0.68) 








To extend the recombination mechanism to the LHC, we need criteria to extrapolate the parameters of 
the model. In absence of data, we compare with the theoretical results obtained with the hydrodynamical 
approach in [12]. We find that, by rescaling v± and tAj_ to the values reported in the second line of 
Table dl a reasonable agreement is obtained for slope and normalization of the predicted yields for vr, p 
and K. We also take unit fugacities for central collisions, although this has a negligible effect in the final 
result. In [26j different values for the extrapolation of the value of v± from RHIC to LHC have been 
found to affect the final result for the suppression. Here we take the conservative approach of fitting 
hydrodynamical results. 

For the calculation of Raa, we model the soft spectrum in the p+p collision by the same Eqs. ([8j)- (fT0|) 
fitted to proton production data from Ref. [6j. This leads to values reported in the Table. The indicated 
fugacities are consistent with indications from the statistical model of hadron abundances; for strange 
particle production we have multiplied by an additional factor 75 = 0.3. The volume rescaling factor 
~ 9 found does not distinguish between a decrease in r and/or A_i_. Just to put some numbers, taking 
instead Rp = 1.2 fm a value of r ~ 1.6 fm is obtained. Finally, the same flow velocity v± is assumed for 
p+p and A + yl collisions at a given energy. These details do not affect much our results since within this 
setup the spectra are, anyway, dominated by the hard component in the pj-range under study. We show 
in Fig. [T] a comparison with experimental data. The agreement for both pion and proton production is 
very good, while it turns out to be worse between our K~ and the data for proton-proton. This can be 
traced back to the behavior of the fragmentation function in the hard part of the spectrum. We did not 
attempt to improve on this, since we prefer to maintain all the parameters fixed in [8j. 

Elliptic flow. It is customary to expand the azimuthal dependence in the function fa, Eq. ([T]) according 
to: 

fa = l + 2v2{p±)cOs{2^a) + -.. (14) 
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Figure 1: The comparison of the spectra of produced tt"*", K and p + p in proton- proton collisions at -^s = 200 GeV with 
STAR experimental data from [B] 



We see at once from Eq. ([7]) that the hadron azhnuthal distribution is: 

/j^ = 1 + 2^2 (P±)cos(2$) + ... (15) 

with: 

V2{Pl_) ^ 2V2{^) (16) 

for small f 2 . For a hadron made of N constituents we recover the well-known scaling law [25] : 

1f2(P^) = V2{^). (17) 

The observation of elliptic flow for the /q and the determination of the value of N by which the 
scaling rule (fT7j) is obeyed would provide a very reliable determination of its quark structure, albeit a 
very difficult one, given the smallness of the expected value of (V2)/o, see ref jTE] for a similar proposal 
for the hypothetical pentaquark. 

Parton fragmentation. For p± large enough, particle production is dominated by hard processes 
computed by: 

a Q 

where Da^h{z,Q'^) is the fragmentation function of a parton a into a hadron h, evaluated at = Pj. 
The corresponding perturbative spectrum of parton a: 

dNr' C 



= K- —5 (19) 

is taken from LO pQCD calculations for central collisions from [27], where the parameters C, i?, and (3 
can be found. A constant K factor of 1.5 is included [8j. 

The suppression of particles with p± > 6 GeV is described by the energy loss of their parent partons 
due to radiative medium- induced processes [3]. We adopt the simplified description in j8], whereby the 
energy loss is implemented by a shift in the perturbative spectrum by the amount: 

Ap^(6,]5^) = e(6)v^^ (20) 
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where: 



e{b) = eo 



l-exp(7^-2) 



1 - exp(-2) 



4 



R\-\ +Ra-\ 



(21) 
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This gives a common suppression for all hadron species, in agreement with experimental data at large 
p_\_ and with the results from [28] . The corresponding value of eo for the LHC is also taken to satisfy the 
results from [28] with a multiplicity enhancement from RHIC to LHC as given from [29j . Notice that in the 
full jet quenching calculation, geometric effects - surface dominated emission - are very important |28j . 
Here we just hide all these effects into the single parameter, eo, which effectively describes the calculated 
suppression. The surface emission makes the extrapolation from RHIC to LHC rather insensitive to the 
actual value of the jet quenching parameter, q, in the approach of [28]. For this reason, we do not expect 
our main results to change due to this uncertainty. 

For the fragmentation functions Da^h{z, Q"^) we use KKP [30j for all the particles except A + A which 
is taken from DSV [31]. These parameterizations are obtained from fits to experimental data and do not 
exist for the H or the /o as a four-quark meson, which will be discussed in the next Section. 

To normalize peripheral or proton-proton collisions to the central ones, we shall need the number of 
binary nucleon-nucleon collisions in Au-|-Au and Pb-|-Pb. The number of collisions for the first case is 
taken from [8]. For the LHC, we have rescaled Ncoii with the inelastic cross section |12j . 

The adopted values for eo and N^oii are given in Table [TJ 

In summary, the total spectrum is just the sum of the soft part described by Eqs. (f8]l- (fT0]l and the 
hard part given by Eqs. (jl8p -(|19p. The relative importance between both contributions depend on the 
quark structure of the hadron and the fragmentation functions. In this sense, the recombination spectrum 
for a 4-quark state is broader than that for a normal meson. We will see in the next section that the 
fragmentation is also more suppressed in the first case than in the second. Both effects make the nuclear 
modification factors Raa and Rcp to be enhanced at intermediate pt for particles with larger number of 
constituent quarks. 

3 Fragmentation functions 

The main uncertainty in the model is the construction of a fragmentation function for /o as a 4-quark 
state. 

The question of how a hard process fragments into the final hadrons is not fully solved theoretically. 
In the approximation of independent fragmentation (jlSp . fits to experimental data by means of GLAP 
evolution equations have been provided. In these fits, experimental data from e^e~ hadrons, in par- 
ticular from hadronic decays are used. One additional difficulty comes from the kinematic constraints 
which, in these experiments, reduce the sensitivity to the fragmentation functions at large values of z. 
On the other hand, for hadronic collisions the relevant values are in fact 2 ~ 0.7 [32] . due to the steeply 
falling perturbative spectrum ()19p . At present, the description of experimental data in p + p collisions at 
RHIC is rather good for mesons and reasonable for baryons. 

We do not attempt a full study of different fragmentation functions and their evolution. Instead, we 
make use of the known fragmentation functions for vr, K, p [30] and A [31] and adopt simple rules to 
extrapolate to the yet unknown case of /o as a 4-quark state, comparing when possible with the scanty 
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available data. We consider the case of H +5"*" as well, for which preliminary data are available from 
RHIC. 

The LEP Collaborations have measured, see e.g. the Z*^ — > /g fragmentation function: 

J_^^ 2Tg j^h^ Ml). (22) 

Thadr dz ^ 3r, + 2r„ " ^ ' ^ ^ 

q=u,a,s 

The fragmentation function for q— > /o (g = u, d) at ~ (5 GeV)^ has been measured in f^-nucleus 
scattering by NOMAD [3^ . 

We make the assumption that the fragmentation functions at large z are dominated by the quark 
content of the particle. In addition, we introduce a suppression factor for each additional strange quark 
with respect to A, which we estimate to be ~ 0.5 from K/t^ fragmentation ratio of KKP ^30j. This leads 
to: 

DT^^\z,Q^)^Dt+\z,Q^), 

since H~ contains d and s but no u valence quarks. 

In the case of the /o as a 4-quark state, we impose an additional suppression factor (1 — z)^'^ . 
Intuitively, this term corresponds to the suppression for producing and additional diquark with z — > 1. 
This suppression factor is present in the numerical comparison between p and vr in the KKP fragmentation 
functions [30], almost independent of the value, and has also been proposed in terms of Regge theory 
[35j . For the gluon— > /o fragmentation function we take a larger suppression factor, 0.1(1 — z)^, although 
this gluon contribution is almost irrelevant. In conclusion we take {q = u, d, s): 

D/«(^«)(z,g2)^ 0.5(1-^)1-5 ^^""^^^^'^'^ 

Z?/o(^'')(z,g2) ^ 0.1(1 - zf ^^^^^(^'Q') (23) 
Finally, for the /o as a ss state, we scale D^'^ from [30j by: 

D{t'\z. Q') - 0.5 X D^Jiz, Q2) . (24) 

We show in Fig. [2]the — > /o(4q) fragmentation (j23p at = GeV^ (blue curve) compared with 
OPAL data [33j (stars). In the same figure, we show zD^°^^''\z) at = 25 GeV^ (red curve) compared 
to NOMAD data (boxes). Within errors, the drop we have assumed for z — > 1 is not incompatible with 
present data. 

Summarizing, the results - to be presented in the next section - are sensitive to the behavior of 
the fragmentation functions at large values of z ~ 0.7 [32j, which control, in particular, the relative 
normalization of the hard and soft spectra for each particle species. Our strategy was to construct the /o 
fragmentation functions from the known functions by applying the rules given above. We have checked 
that the results are still compatible with the small amount of available experimental data. We have a 
check on the quality of this procedure in the case of the H, where no fragmentation functions are available 
and the description of the data, applying the same rules as for the /o, is good - see Fig. O 
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Figure 2: the Z° /o(4q) fragmentation at = lO" GeV^ (blue curve) compared with OPAL data [33] (stars). In 
the same figure, we show zD{°'^^''\z) at = 25 GeV^ (red curve) compared to NOMAD data (boxes). The drop we have 
assumed for z — ^ f is not incompatible with present data, within errors. 



4 Results 



We present in Figs. I3I4I and Fig. [5]our results for the nuclear modification factors Rcp and Raa at RHIC. 
The nuclear modification factors are defined as: 



Rcp 
Raa 



A^coll(^) (fN^^+Anih = 0)/dPl 
N,^ll{b = 0)dmAn+An{b)/dPl 
d''NAu+An{b=0)/dPl 



N,oii{b = 0)dmp+p/dP]_ 
We take b = 12 for peripheral collisions at RHIC. 




(25) 
(26) 



Figure 3: Rcp, Eq. ([25}, for n" and p+p in Au+Au 
at RHIC. Peripheral collisions at 6 = 12 fm. Data 
from l36l. 




Figure 4: Red curves: Rcp for Ks and A + A, 
Au+Au at RHIC. Peripheral collisions at 6 = 12 fm. 
Blue curves: predictions for /o = ss (dashed) or /o = 
[qs][qs] (full). Data from [32]. 



For vr, protons, K and A we reproduce the results of Ref. [SJ. The large difference between Rcp 
and Raa for the strange baryons is, in this approach, identified to be due to the strangeness suppression 
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factor in p + p, as it was also found for the total yields at SPS energies [14^ I15j . and in agreement with 
PI- 




Raa 



flit 



2.5 



1.5 



0.5 



10 




Figure 5: Starting from below: Raa, Eq. (f26|l . for vr, Ks, /o(980) as a ss state, p + p, A + A, /o(980) as a 4-quark state 
and + E'^ . Data from Ref. [38) . 

Figs. [4] and [5] illustrate clearly our main result, namely that there is a large difference between the 
suppression for the /o(980) when assumed to be a four-quark state or a quark-antiquark meson. For 
both Rcp and Raa, the nuclear modification ratio for the 4-quark state is in the range of the strange 
baryon ratio, while it is closer to the K5 signal for a ss state. The difference is particularly large in the 
experimentally accessible range oi p± < Q GeV. 

Taking our calculations at face value, the effect is larger in the case of Raa, but the uncertainties are 
also larger, due to the implementation of the unknown fragmentation functions for these mesons, which 
in the case of Rcp cancel to some extent. Nonetheless, even if the effect could be overemphasized in the 
present formalism, we expect Raa for the four-quark state to be in the range of the corresponding values 
for A -|- A and H~ -|- H^. This gives still a large leeway to distinguish the valence quark structure of the 
/o(980). 





Figure 6: Rcp for LHC using the same notation as 
in Fig. E] 



Figure 7: Raa for LHC using the same notation as 
in Fig. H 



To extend the formalism to the case of the LHC we extrapolate the parameters as discussed in 
Section [21 both for the hard part of the spectrum, given by fragmentation, and for the soft part. For 
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peripheral collisions we take 6 = 14 fm. The corresponding predictions with this set of extrapolating 
parameters are presented in Figs. [6] and [71 

The Rcp and Raa signals are considerably reduced, presumably due to the large increase of gluon 
and sea-quark densities in the proton, which entails an increase in the fragmentation component. This 
conclusion, however, is strongly sensitive to the assumed value of the recombination volume, tA±, of the 
fugacities at LHC as well as to the values of v± ^6]. Cross section data are needed for a more reliable 
assessment of the parameters of the model, as is the case at RHIC. Nonetheless, the correlation of the 
^CP/AA ratios of a 4-quark /o with the baryon ratios is still borne out by our calculation. 

5 Discussion 

We have presented experimentally observable quantities, which in the region of intermediate transverse 
momentum in heavy ion collisions are sensitive to the constituent quark content of the resonances. For 
this purpose, we use a simple but successful model based on a hadronization by recombination of the 
quarks in the fireball produced after a heavy ion collision [8] . We extend this model to the multi-strange 
and multi-quark hadrons, finding a good description of both Raa and Rcp with the usual assumptions 
on strangeness suppression. We extrapolate the model to the LHC by fitting to the available results in 
more refined calculations of both the soft [12] and the hard contributions to the spectrum. Uncertainties 
however remain in the values of the parameters at the LHC, which will be possible to eliminate when 
cross section data will become available. 

The /o(980) resonance is, at present, the best possible case for these type of studies as (i) the signals 
are large enough for an unambiguous result to be obtained and (ii) the resonance is abundantly enough 
produced [18] to afford a measurement of the nuclear modification ratios. 




P- (GeV) 

Figure 8: Rcp for LHC. We compare the two hypotheses, 2(7-dashed and 4(j-solid, for the /o(980) and /o(1370). The four 
quark interpretation of /o(1370) is q being a hght non-strange quark [21] 

Other possibilities concern the /o(1370) and /o(1700) mesons, for which reasons to suspect a 4-meson 
interpretation have been recently advanced [21], alternative to the conventional interpretation as P-wave 
qq states [20] . We do not have any information on fragmentation functions into these mesons analogous 
to those provided by OPAL and NOMAD, see Fig. [5J To give an idea anyway, we simply scale with 
the strange quark content, taking for /o(1700) the same as in (j23p and multiplying by 1/(0.5)^ for 
/o(1370) [21]. We report in Fig. [5] the resulting values of Rcp at RHIC of the /o(1370) according to 
the two assumptions on composition, compared to the analogous curves for /o(980). The predictions for 
/o(1700) essentially coincide with those for /o(980). 
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In the same mass region, the /o(1500) scalar meson has been long suspected [20] to be the lightest 
glueball predicted by QCD lattice calculations (see e.g. [39] )• We may advance two possibilities. One 
is that recombination of constituent gluons exists. Since the lightest glueball should be made by two 
constituent gluons, the /o(1500) would then follow a nuclear modification pattern close to that of con- 
ventional mesons. Alternatively, glue recombination could be suppressed and the /o(1500) signal arise 
from its mixing with the /o(1370) and /o(1700) mesons only. The /o(1500) would exhibit, in this case, 
an anomalous suppression in the coalescence region. 

Experimental investigations at RHIC and LHC could shed light on the above issues and produce a 
decisive advance in our understanding of meson spectroscopy. 
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